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ABSTRACT 

X-ray telescopes having a relatively wide field-of-view and spatial resolution vs. polar off-axis angle curves much 
flatter than the parabolic dependence characteristic of Wolter I designs are of great interest for surveys of the 
X-ray sky and potentially for study of the Sun’s X-ray emission. We discuss the various considerations affecting 
the design of such telescopes, including the possible use of polynomial mirror surface prescriptions, a method of 
optimizing the polynomial coefficients, scaling laws for mirror segment length vs. intersection radius, the loss of 
on-axis spatial resolution, and the positioning of focal plane detectors. 
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1. INTRODUCTION 

X-ray optics involving polynomial extensions of Wolter I mirror prescriptions are of astronomical interest in the 
context of wide-field X-ray surveys, the goal of the proposed Wide-Field X-ray Telescope (WFXT) mission, and 
potentially also for study of the Sun’s X-ray disk. In this paper, we explore methods for designing, optimizing, 
and evaluating such optics. 

Mirror presciptions and ray trace methods are reviewed in §2 — §7. Expressions for evaluating spatial reso- 
lution and effective area using Monte-Carlo simulations are presented in §8 and §9. The approximate analytic 
expressions for Wolter I optics given in VanSpeybroeck and Chase 1 are reviewed and extended to a broader 
domain in §10. Effects due to figure errors and surface roughness are briefly discussed in §11. Our methods for 
designing and optimizing telescopes are described in §12, including the topics of figures of merit §12.1, design 
considerations §12.2, single shell optimization §12.3, and nested configuations §12.4. A few closing remarks and 
plans for the future are presented in §13. 

2. MIRROR PRESCRIPTIONS 

We assume the lengths and nominal grazing angles for the primary (P) and secondary (H) mirror segments are 
the same. We place the origin of the optical (z) axis at the intersection plane between the two mirror segments, 
with z increasing toward the front of the primary and decreasing toward the back of the secondary. (This points 
the z-axis in the direction opposite to that in Burrows, Burg, and Giacconi, 2 henceforth BBG).) The mirror 
surface prescriptions for the conical approximation are: 


rl tP {z P ) = rl [ 1 + tan (a 0 ) {zp/r 0 )f 

(1) 

t I,h( z h) = rg [ 1 + tan (3a 0 ) (z H /r 0 )] 2 

The mirror surface prescriptions for Wolter I and polynomial optics are: 
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( 2 ) 


r l,p( z p) = r l [ 1 + 2 tan(ao) ( z P /r 0 ) + 6a 2 ( z P /r 0 ) 2 + (5a 3 {z P /r 0 ) 3 ] 

r l,H( z p) = r o [ 1 + 2 tan(3ao) (zp/ro) + tan 2 (3ao) /i(ao) (.zp/ro) 2 

+ Sb 2 ( zp/r 0 ) 2 + &> 3 {zp/r Q ) 3 ] 

fe(ao) = 1 - 1/(1 + 2cos(2a 0 )) 2 , 

where 

«o = (i) tan” 1 Cj), (3) 

ro is the radius of the segments at their intersection plane ( z P = zp = 0), and / is the the distance from 
the intersection plane to the on-axis focus for Wolter I optics. Thus the on-axis focus for Wolter I optics is at 
z = —f. We denote the mirror segment length by l. The mirror shell length is 2£ Wolter I optics require 
Sa 2 = < 5(23 = Sb 2 = Sb 3 = 0 . 

3. APERTURE FOR INCIDENT RAYS 

For a given shell s, the maximum radius of the aperture, r moXiS , is independent of off-axis polar angle, 6, and 
position azimuthal angle, <j>. The value of r rnax , s is determined from the equation for the surface of the primary 
with 2 set equal to L The minimum radius of the aperture, r m j„ iS (0, 0) , is a function of 6 and <j>. Straightforward 
geometry leads to the expression 

r m in,s(0, 0) = -t tan0 cos^> + \Jrl s - i 2 tan 2 8 sin 2 0, (4) 

where ro, s is the intersection radius for shell s. Amusingly, although the shape and centroid of the incident 
aperture varies with 8, the total geometric area of the incident aperture is independent of 9, and is given by 

Ageom,s 2T ( r moa , )S ^0,s )' (^) 

These expressions neglect possible shadowing by the next inner shell, a pre-collimator, or any support structures. 

4. RAY INTERSECTION WITH OPTIC 

The equations of ray tracing X-ray optics are simple in principle. Given a starting position (xq, y 0 , Zq) and 
normalized wavevector ko = (k x .o, k y .o, fc z ,o), the intersection (aq, yi, zi) with an optic is determined from: 

Xi = Xq -{- kx,Q t 

y\ = Vo + fey, 0 t (6) 

Zi = Z 0 + k z> 0 t. 

The parametric variable t is determined from 


( Xq + fe x , 0 t ) 2 + ( y 0 + fey, 0 t ) 2 = Tq{zq + fe z ,o t ). 


( 7 ) 



Here ro(z ) is the radius of the optic as a function of distance along the optical axis. For simple ideal surfaces with 
no figure errors, ro(z), Eq. (7) can be solved analytically for t, and hence for (xi,yi,zi). For more complicated 
surfaces or in the presence of surface errors, we proceed by first finding the intersection (xio, t/io, .zio) with a 
simple surface, ro,o(z), approximating the actual one, r(z) = ro,o(z) + A r(z), and then correcting the solution 
using a perturbation method. The 1st order corrections to the intersection point are given by: 


X\ — Zio “I” &x,0 St 

2/i = 2/io + k y fi St (8) 

zi = z w + k Zt o St, 


where St is given by 


St = f 0 , ,(*..) Mzin) (9) 

fcx ,0 2? 10 H“ fcy , 0 2/10 fcz , o ro,o(zio) 

dz U=*10 

In Eqs. (8) & (9), we’ve neglected radial perturbations dependent on x and y. This is justified because contribu- 
tions from axial slope perturbations dominate those from ( x , y, z) perturbations in intersection points or from 
azimuthal slope errors to the PSF. This method can also be used to deal with low-frequency figure errors. 

5. SPECULAR REFLECTION FROM A SURFACE 

Neglecting azimuthal variations as is appropriate for ideal surfaces, unit normals ni = (n X) i, n Vi i, n z>1 ) to the 
surface O at the intersection point (*i, yi, zi) are given by 



where 



For example, for a polynomial primary surface: 


dr P}P _ tan(q 0 ) + <kt 2 (z/r 0 ) + ( 3/2)fia 3 (z/r 0 ) 2 
dz y/l + 2 (z/r 0 ) tan (a 0 ) + Sa^z/ra ) 2 + 5a 3 (z/r 0 ) 3 


(10) 


( 11 ) 


(12) 


The direction vector k\ = (k x .i, k y .i, k z .i) of the ray specularly reflected from the surface is given by 



k\ = (hi x ho) x hi — hi (ko ■ hi) = ko ( hi ■ hi ) — 2 hi ( ko ■ hi ) 


fc*,i = k x> o (n^j + - n^ j) - 2 n X]i (k yfi n yA + k zfi n z>1 ) 

k y , i = fcy,o (n^,i + n*,i - n^j) - 2 n y> i (fe s>0 n x ,i + k zfi n z , i) 


(13) 


kz,l — k Z: 0 (^ Xj i “1“ Tl y ^ i 2 7l 2) l (k x , 0 «*,! "f“ &y,0 Tly, l)* 

The right-hand-side of the top relation in Eqs. (13) follows from the vector identity (A x B) x C = B(A ■ C) — 
A(B ■ C). The vector hi will be normalized to unity if ko and hi are, as the reader can demonstrate for himself. 

If included directly in the ray trace, scattering due to dust or surface roughness introduces perturbations to 
the specularly reflected direction vector. 


6. PROPAGATION TO THE FOCAL PLANE 

Given the intersection position ( x 2 , yi, z?) with the secondary and the normalized wavevector (k x .2, k Vt 2, and 
k z o) reflected from the secondary, the ray is propagated to the focal surface plane 2 = — / + 5 z by solving 
— / + Sz = Z2 + k z _2 t for t and substituting in Xf = x 2 + k Xi 2 t and y/ = 3/2 + k Vi 2 t to find Xf and yf. 


Xf = X2 + “ (-/ -1 - Sz - Z2) 
Kz,2 

1 if = yi + ^ (-/ + Sz - z 2 ). 

Kz,2 


(14) 


From Eqs. (14), and the quantities X2, 2/2, 22, 2/^,2 and 2/^,2 we can construct the rms blur radius, 

&rms (fit Sz), as described in § 8.1. 


7. GRAZE ANGLES 

The graze angles, ai and 0:2, for the ray at the primary and secondary, respectively, are: 



Values of the graze angles are needed for applications requiring the calculation of reflectivities for the weighting 
of individual rays. 



8. MEASURES OF SPATIAL RESOLUTION 


8.1 Calculation of cr rms from a ray bundle 

RMS dispersion radii for PSFs derived from a Monte-Carlo simulation are determined as follows: 


<Q> = 



(16) 


—£> I A 

° x.rms ' ^ x.rms 


N 


N- 


y(< X > - <X> 2 ) 


N 


a,, 


N 


— (<y 2 > - < y > 2 ), 


(17) 


where N is the number of doubly reflected rays reaching the focal plane. We define Xi = Xfj + Sz(k Xt i/k z .i) and 
Ui = Df,i + Sz(k y .i/k z .i), where Xf.i is the ray ^-position in the fiat focal plane through the on-axis focus, k x ,i 
the ^-component of the doubly reflected ray direction vector, Xi the ray ^-component at a height Sz above that 
plane, and similar definitions for the '(/-components. Then Eqs. (17) become: 


^rmat^) ff i,m»(^) ^ y.rms (^^0 > 


(18) 


— (a l) ( < " x > r7ns "f" ^ b x ,rms ®i,raa ) 


< Xf > - < Xf > 


(fe) 


< */ l IT J > - < */ >< l Y I > 


. k x \ ( k x 

< 'fc > - < U l > 


(19) 


® y,rms (^0 
Cy,rms 


1 ^ (py> rrns ^ by^rms fiz “ 1 “ Q'y,rms ) 


< Vf > - < Vf > 2 


vy^rms 


(fe) 


< Vf hr > - < Vf >< hr > 


ky 




( 20 ) 


For a flat detector with face perpendicular to the optical axis, the optimal focal height, 6z m i n , determined by 
minimizing the RMS blur a rms , then is: 



Sz • — — 


bx,rms by,rms 

( Q> x,rms Q , y,rms ) 


Substituting Eq. (21) into Eqs. (19) — (20), and assuming N/(N — 1) ~ 1, leads to: 


— Cx,rms ^y,rmsi 


( 21 ) 


( 22 ) 


Jl (x~ \ _ t(\\ (fix,rms + ^y,rms ) 2 / 0 „s 

\n X) rms i Uy^rms) 

These equations apply independent of the mirror prescription. 

For Wolter I optics, 8z m i n is nearly a parabola in the polar off-axis angle 0 , and the best focal surface is 
usually not a flat plane. This leads to schemes for arranging focal plane detectors to follow the best focal surface 
as closely as physically possible. 

8.2 Measures of encircled energy 

Encircled energy ( EE(r )) vs. r is just the fraction of rays falling inside a radius r, and obviously varies from 
0 at r = 0 to 1 as r -too. Common measures of spatial resolution are rso = HPD/ 2, the radius inside which 
50% of the rays fall, and rso = £>8o/2, the radius inside which 80% of the rays fall. Here HPD is the PSF 
half-power diameter. For a two-dimensional gaussian PSF with a 1/r cusp, one finds HPD = 1.35 a rms and 
.DgO — 2.56 (Tj 'ms- 


9. EFFECTIVE AREA A E ff{9) AT E = 0 
In the zero energy limit, where reflectivities are all unity, the effective area is defined by 

Ae ff {E = 0, 0) = A geom (0) x n 2 (0)/ n inc (0), (24) 

where A geom (0) is the geometrical area of the incident aperture (whether for a single shell or a set of nested 
shells) at polar off-axis angle 0 (see § 3), n 2 (0 ) the number of doubly reflected rays through the optic, and rii nc (0) 
the number incident on the aperture. 

10. APPROXIMATE ANALYTIC EXPRESSIONS FOR WOLTER I OPTICS 

We now consider Wolter I optics, for which Sa 2 = = 5b 2 = db-j = 0 in Eqs. (2). 

10.1 Single shell 

VanSpeybroeck & Chase 1 provide a number of useful approximate analytic expressions characterizing the per- 
formance of Wolter 1 optics. For the RMS blur radius, u rms .o{0), as a function of off-axis angle, 0, in a flat plane 
perpendicular to the optical axis at the on-axis focus ( 5z = 0), their expression is 

<^rms,o(9,ao, fj,o) — IM) (2/4) (i/f ) (tan 2 0 / tana o) + 4 tan0 tan 2 ao. (25) 

In Eq. (25), ao = (1/4) tan -1 ( ro/f ), the graze angle at the optic radius of the primary-secondary intersection, 
/ the focal length, l the mirror segment length (so that 2 t is the shell length), and jig ~ 1 an adjustment 
factor. The second term on the RHS of this equation is due to coma. Since coma grows linearly with 0, while 
the spherical aberration term grows as 0 2 , coma is only important at the smallest values of 0. In the flat on-axis 
focal plane, the cross-over occurs at tan# ~ (8 / / £) tan 3 ao- 

The VanSpeybroeck & Chase 1 expression for the best focal position, Szb es t(9), above the on-axis focal plane 
as a function of off-axis angle, 0, is 



Szbest(0,a o ,fJ.Sz) — Viz (l/ 8 ) ^ (tan0 / tana 0 ) 2 , (26) 

where f. i$ z ~ 1 is another adjustment factor. On the best focal surface given by Eq. (26), VanSpeyBroeck & 
Chase 1 find that the RMS blur radius, (j rrns (8 ) , neglecting coma, is 

O rms,best{0,a o ,fH ) ~ [h (1/4) (e/f) (tan 2 0 / tana 0 ). (27) 

Comparing Eqs. (25) and (27), we see that on the best focal surface the first term on the RHS, due to spherical 
abberations, is half its value in the on-axis focal plane. 

In the low energy limit E — ► 0, VanSpeybroeck & Chase (1972) find that the effective area as a function of 

8 is 


A-ef f[8 ,CJLQ, 1 A geom [1 f^A,2 (2/3) [8 f a 0 )] , (28) 

where A geom is the on-axis geometric area, and ha.i ~ 1 and fiA , 2 ~ 1 are adjustment factors. 

10.2 Single shell generalization 

The results of §8.1 show that we can write the RMS blur radius, in linear dimensions (such as cm), as a function 
of off-axis angle 8 and height Sz above the flat on-axis focal plane as 

a rms(°’ Sz ) - C S + 2 b s Sz + a s 8z 2 . (29) 


Here 


Cg = < Xf + y 2 f > - < Xf > 2 - < y f > 2 



According to the results from §10.1, we can write 

Cg(8,ot o,/i c ) — [ A*c (2/4) l (tan^O / tana 0 ) +4 / tan8 tan 2 a 0 ] 2 , (31) 

where p c rsj \ . 

Monte-Carlo simulations show that a very good approximation for a s is given by 


a s («o,/-ta) — (4 /i Q tana 0 ) 2 j 


(32) 


where /i a rsj \ . 

The condition that the first term in Eq. (25) is cut in half on the best focal surface now gives us a way to 
determine b s . Neglecting the coma term in the expression for c s , the result is 


b s (6,fi c ,fi a ) ~ - n a n c \/3 l tan 2 8. 


(33) 



So far this expression has been followed by the results of the Monte-Carlo simulations for Wolter I optics that 
we have done. Then we have 


o^ms{Sz,0,ao >a ,fic,tia) ~ c s (6,a 0 ,fi c ) + 2 b s (6, fi c , fi a ) Sz + a s (a 0 , fi a ) Sz 2 

(34) 

% ~ b s(@ , /^ci l^a) 

vZmin — 7 \~ ■ 

0 's [ OLq , ft a ) 


Substitution of a s (ao,fi a ) and b s (6, fi c , /i a ) in the last expression above reproduces the dependences in the 
expression for the best focal surface, Eq. (26). 

10.3 Nested shells 

Eqs. (29) — (34) are for a single Wolter I shell. For a telescope consisting of nested shells, we might consider 
summing expressions for each shell, weighted by effective area: 


N S 

Aeff,nested{@i /^4i = ^ ^ ^-ef f i ^0,s> M-A,2)> ( 35 ) 

S=1 


&rms,Q, nested {@1 ftci ftA,li ft A, 2 ) 


finest, nested {@1 ftci V ai ft A, li ft A, 2 ) 


&rms, 


nested {@1 ftci ftai ftA,h ft A, 2 ) 


( 

\-A-eff, nested ($> ft A A , ft A ,2) 

N S 

V. Aeff (6,a 0>a ,llA,l,HA, 2 ) Grms (0, 0, O!0 jS , ftci fto) 
8=1 

( I 

\A e ff, nested (fit ft A A > ft A , 2 ) 

Ns 

^ ) Aeff ($> C*0,s > ftA,h ftA, 2 ) bZb es t( 9 , C*0,s> ftci fta) 
s=l 




( 1 ) x 

\A e ff } nested (^5 ft A, 1 , ft A, 2 ) / 

Ns 

^ ^ Aeff (@: O:o,s ; ftA,li ftA,2) (?rms(Sz,Q ,Cfco,S) Me? Mo) 
3=1 


(36) 


11. EFFECT OF FIGURE ERRORS AND SURFACE ROUGHNESS 

All equations have been for ideal figures. There will of course be figure errors and surface roughness resulting 
from the manufacturing process. Once the optics are manufactured and their figures measured, it is possible to 
evaluate the effects of sufficiently low amplitude and low frequency figure errors using the perturbation equation 
given in §4. It is also possible to evaluate the effects of surface roughness by convolving the predicted scattering 
wings with the predicted geometric core (see, e.g., O’Dell et al.(1993) 3 ). A rough estimate of the resulting 
dispersion follows from the root-summed-square of the contribution expected from manufacturing errors and the 
contribution expected from ideal figures. We have already mentioned that positioning focal plane detectors to 
follow as closely as physically possible the best focal surface leads to the best performance. 



12. WIDE-FIELD X-RAY OPTICS WITH POLYNOMIAL COEFFICIENTS 


In 1992 BBG 2 showed how by adding higher order polynomial terms to the Wolter I prescriptions, and hence 
giving up some on-axis spatial resolution, one can obtain prescriptions for reflecting surfaces that provide im- 
proved spatial resolution over a wide field-of-view (say ~ 30 — 45 arcmin). Such so-called polynomial optics, 
Eqs. (2) with non-zero 5a,2, 60,3, 8b 2, and <5 63, would be particularly useful for moderately deep to deep surveys, 
to be carried out by observatories such as for the proposed Wide-Field X-ray Telescope (WFXT) mission, and 
for solar x-ray observations. 

12.1 Figures of merit 

Many of the methods outlined above for Wolter I optics carry over to the case of polynomial optics. The most 
important new ingredient is how to choose the polynomial coefficients 8a 2, Sa 3) £62, and Sb 3 so as to optimize the 
telescope performance for the chosen task. In the case of a survey mission, for example, in principle one 
needs to create many simulated fields-of-view, provide accurate background, optical and detector 
performance models, and chose the polynomial coefficients and other mission parameters, subject 
to spacecraft and launch constraints, that unambiguously find the most sources. In practise, one 
often designs the optics to optimize some suitably chosen merit function, 2, 4-8 and then test its performance using 
simulated fields-of-view. Methods for opimizing parameters can get complex, one example being the Markov- 
Chain-Monte-Carlo method. 7,8 

One simple merit function stressing the minimization of the spatial resolution over some specified field-of-view, 
Ofov, is 


M(8z, 6 fov , 8a2 , Sa3 , Sb2 , Sb3 ) 


Jq FOV dO sin 0 w{6 ) a? rns {8z(0),0,8a 2 ,8a 3 ,8b2,6b 3 ) 
Jq FOV dOsinO w{6) 


(37) 


The weighting factor, w(6), is often set unity, but may also be set to the effective area, A e ff(6). Since maximizing 
A e f f (0) is also important for greatest sensitivity, one may want to consider the weighting factor w(0) = 1 /A e f j (0) . 
However, normally the effective area is constrained mostly by mass considerations and spacecraft dimensions. 
Eq. (37) provides a spatial resolution figure of merit for the full field-of-view. It is important to recognize that 
polynomial optics with ideal surfaces necessarily results in the loss of on-axis spatial resolution. 
Of course, manufacturing errors may make this point irrevelant. 


12.2 Design considerations 

Since X-ray optics necessarily involve grazing incidence, effective area is built up by nesting many shells in each 
telescope module. This complicates the optimization procedure because: 


1. The polynomial coefficients for each shell will be different, depending on the shell intersection radius, ro, s , 
and segment length, £ s . Here s denotes the shell number; our convention is to number the outermost shell 
as s = 1 with s increasing inward. Most telescope designs for wide-field surveys contain many shells (as 
many as 100) , so this fact greatly complicates the optimization procedure. In practise one often optimizes 
the polynomial coefficients for each shell separately rather than all at once. 

2. The best focal surfaces for the different shells will be different depending on the shell segment length, l s . 
Thus the optimization procedure should include the possibility that the segment lengths of the various 
shells are not all the same. The optimization procedure should also include placement of the focal plane 
detectors. 

3. Constraints due to mass limits, spacecraft dimensions, and optic manufacturing limits normally lead to an 
envelope for the telescope optics prescribing values for the maxima and minima intersection radii and shell 
Segment lengths 7*0 : rnax and 7* O.-min , and £m(ixi and ■ 



For Wolter I optics, we see from Eq. ( 26 ) and the fact ccq = ( 1 / 4 ) tan 1 (ro//) ro/ 4 / (see Eq. ( 3 )) that 


Szbest( 0 ,a o ,fisz) - Viz ( 1 / 8 ) l (tan# / tanao) 2 ~ fj, Sz 2 / ( 4 ) # 2 ( 38 ) 

r o 

Eq. ( 38 ) shows that Wolter I mirror shells satisfying the scaling relation 

ro >s 

TO, max 

are expected to have approximately the same best focal surfaces. The difficulty with this scaling relationship 
is that for closely packed shells, necessary to build up effective area and to reduce or eliminate ghost rays 
due to single reflections, the limit on shell length is normally reached well before the intersection radius ro, s 
approaches the value ro, m in, thus limiting the effective area. In addition, of course, introduction of the higher 
order polynomial terms in the surface equations alter the shape of the best focal surface. 

This suggests one should investigate scaling laws of the form 

(-max > t.=tmax (-^-)" > (-min- ( 40 ) 

T*Q {max 

If the available space is to be filled up with mirror shells, one has 


'-max — 


> L = ir 


) 2 > 4 


( 39 ) 


^ log {(max / (min) ( 41 ) 

log (r 0 ,max /r 0 ,min ) 

12.3 Symbolic optimization of a single shell using Mathematica © 

We use the symbolic mathematics system Mathematica® 9 to get an expression for M ( 6z , 0 F ov, $ a 2> &I3, &>2, <^>3) 
for given numerical values of ro , / (and therefore c*o) and (. This is done at each stage of the ray trace (intersection 
with mirror surface, computing normals and reflected direction vectors, etc.) by expanding the appropriate 
quantity to second order in da 2, &13, 6b 2, and < 5 63, including cross-terms. For example, the ^-component of the 
normal vector from the secondary (H) surface is computed in the form: 


n x ,H = n X H, 0000 + 6a-2 n X H.iooo + 61 13 n X H, 0100 + <^#2 n X H, 0010 + ^63 n x jj ; oooi 


+ 6a 2 n X H , 2000 + d a 3 n xH, 0200 + fib* n X H , 0020 + ^3 n xH, 0002 


+ 6a2 6a3 n X H, 1100 + $ a 2 <5f>2 n X H, 1010 + &&2 6b3 n X H, 1001 


( 42 ) 


+ 6 a 3 6 b 2 u x h.oiio + 6 a 3 (5 63 n X H, 0101 + 6 b 2 <54 ti x h.ooii, 


where each of the 15 n X H,ijki is a number. The normal vectors and direction vectors are kept normalized 
to 1 to second order. Using Mathematica ® , ray traces over a set 0 of 6 values produces expressions for 

(see Eqs. ( 19 ) and ( 20 )) in the form of Eq. ( 42 ). Then we con- 


®x,rms) ^x,rmsi Cx,rmsi & y,rmsi ^y^rmst Cy, 

struct an expression for a xms (6z,8,6a2,6a3,6b2,6b3), which can be integrated over 6 to produce an expres- 
sion for M(6z,9FOV,5a2,6a3,6b2,6b3). Finally the merit function can be minimized by setting derivatives of 
M(6z,9FOVi5a2,6a3,6b2,6b3) to zero and solving for 6z, &2, 6a3, 6b2, and <563. 


Maintaining symbolic expressions as described above is time and memory consuming. Because of these 
limitations, we proceed in three steps. The first step is not a Monte-Carlo simulation at all. Instead, we choose 
rays spaced uniformly, but not randomly, over the full shell aperture. For each value for 8 < Ofov , the subset 
of rays incident on the shell are traced through the optic. 



The second step uses the results from the first in order to determine trial values for 5z, 60 , 2 , 5 a 3, 662, and 663 
for the value of Ofov- 

The third step is a full Monte-Carlo ray trace using these trial values with large numbers of rays per 0 value. 
The only approximation in the full Monte-Carlo ray trace involves calculating the intersections with the mirror 
segments as if they were Wolter I optics, and then correcting the location of the intersections using a pertubation 
method. The effects of this single approximation are insignificant compared to those of, say, axial slope errors. 
This third step provides a full evaluation of the mirror shell performance for the given values of the polynomial 
coefficients. 

12.4 Nested Mirror Shells 

As previously remarked, X-ray telescopes usually involve nested mirror shells. For example the Chandra X-ray 
Observatory has 4 shells and ESA’s XMM-Newton 58. Since the plate scale and best focal surface for each shell 
are different, it is not sufficient to optimize each shell separately. Assuming two coefficients per mirror segment, 
one would need to simultaneously optimize 16 and 232 coefficients for the two telescopes mentioned above. Thus 
far we have succeeded in optimizing polynomial coefficients for 7 mirror shells simultaneously. For telescopes 
with significantly larger number of shells, as currently proposed for WFXT, we have not so far found it possible 
to simultaneously optimize all the polynomial coefficients. 

Instead, for a large number of mirror shells, we proceed as follows: 

1. Given a focal length /, a desired field-of-view Ofov, values for ro. max , ro, m j n , t max , and £ m in, shell thickness 
t, the dependence of £ s on ro )S , and an algorithm for nesting the shells, we produce a Wolter I telescope 
design. 

2. We select a subset of these shells, separated for example by a constant amount in ro, s , for which we 
individually determine the polynomial coefficients &12, &X3, $>2, and 5h^ as described in §12.3. 

3. The coefficients determined in this way form smooth functions of ro, 8 , which we use to interpolate values 
for the polynomial coefficients of the intervening shells. 

4. Now having polynomial coefficients for every shell, we run large numbers of rays in Mathematica© for 
each individual shell. 

5. We use the results of these ray traces to evaluate the performance of the wide-field X-ray telescope we have 
constructed. 

6. One way to improve spatial resolution performance is to shift the individual shells relative to each other 
by a small amount in the axial direction. 2,4,5 

7. A particularly effective way to improve spatial resolution performance is to tilt focal plane detectors so as 
to conform more closely to the best focal surface. This was done with the ACIS CCD detectors on the 
Chandra X-ray Observatory. 

In this way, we are able to construct and fully evaluate a wide-field X-ray telescope consisting of nested polynomial 
optics, as are needed for the WFXT mission. 

13. CONCLUDING REMARKS AND FUTURE WORK 

We have outlined our methods for designing and evaluating wide-field X-ray telescopes for use on X-ray astronomy 
survey missions, with the proposed WFXT mission particularly in mind. Along the way we have extended to a 
broader domain the analytic approximations for Wolter I optics given in VanSpeybroeck and Chase. 1 

Our method for determining the polynomial coefficients for wide-field grazing incidence optics uses the sym- 
bolic mathematics system Mathematica© . 9 This method does not require a Monte-Carlo search of parameter 
space to find values for these coefficients. However, Monte-Carlo simulations are still required to fully evaluate 



performance. Up to the present, we have been able to simultaneously optimize coefficients for up to seven mirror 
shells. Most designs require many more shells and for those cases we still must optimize each shell individually. 
The resulting differences in plate scale and best focal surface can be partially compensated for by varying the 
mirror segement length, £ s , with intersection radius, ro, s , adjusting the relative axial positions of the individ- 
ual shells with respect to each other, 2, 4,5 and tilting the focal plane detectors to most closely conform to the 
telescope’s best focal surface. 

In the future, we intend to: 

1. Explore and evaluate wide-field designs for the WFXT mission. 

2. Investigate designs satisfying the scaling relation Eq. (40) in order to reach the most favorable compromise 
between best spatial resolution over the field-of-view and largest effective area for the available physical 
volume. 

3. Explore incorporating detector tilt into the optimization of polynomial coefficients. 

4. Provide point-spread functions for individual designs in order to enable their evaluation using simulated 
fields-of-view. 
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